Inhomogeneous magnetic phases: a LOFF-like phase in Sr3Ru207 
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Abstract 

The phase diagram of Sr3Ru207 contains a metamagnetic transition that bifurcates to enclose 
an anomalous phase with intriguing properties - a large resistivity with anisotropy that breaks the 
crystal-lattice symmetry. We propose that this is a magnetic analogue of the spatially inhomo- 
geneous superconducting Fulde-Ferrel-Larkin-Ovchinnikov state. We show - through a Ginzburg- 
Landau expansion where the magnetisation transverse to the applied field can become spatially 
inhomogeneous - that a Stoner model with electronic band dispersion can reproduce this phase 
diagram and transport behavior. 
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Fulde and Ferrell [l| and Larkin and Ovchinnikov [21] conjectured that the transition 
between superconducting and insulating behavior, driven by a magnetic field, could occur 
via an intermediate phase with spatially modulated superconducting order. This proposal 
has since been extended to a wide range of settings, from ultracold atomic Fermi gases [3fl 
and exciton insulators m to quark matter and neutron stars I5J. However, experimental 
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confirmation of these predictions is still controversial |fj, ]7|]. In a similar spirit, intermediate 



phases between a Fermi liquid and Wigner crystal [8| have been discussed. We propose an 
inhomogeneous magnetic phase that can be considered a magnetic analogue of the LOFF 
phase. In this case, a change in homogeneous ferromagnetic order occurs via an intermediate 
phase with spatially modulated magnetization. This phase would generate clear experimen- 
tal signatures. Furthermore, we argue on the basis of both new and prev ious experimental 



results that the anomalous phase behaviour observed in Sr 3 Ru20 
explained in this way. 
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FIG. 1: The phase diagram of Sr3Ru207 as inferred from in-plane transport properties. The green 
planes correspond to abrupt changes in resistivity as a function of field. Blue shading indicates 
regions where the in-plane resistivity is anomalously high, becomes highly anisotropic with respect 



to the in-plane component of the field 111 ], and shows an anomalous temperature dependence. The 
phase diagram obtained from magnetic susceptibility [10]] shows the same first order transitions as 
indicated here in green, but lacks the roof. 



The bilayered ruthenate S^Ru^C^ shows a sequence of metamagnetic transitions 12], 
Recent ARPES data has found evidence of van Hove singularities that may drive this meta- 
magnetism 15]. Early studies focussed on a line of metamagnetic critical end-points that 



could be tuned to a quantum critical point by adjusting the magnetic field strength and 
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orientation [13|. Subsequently, ultra-pure samples showed a bifurcation of this line upon ap- 
proaching the putative quantum critical point [si, [lQ] with a second line of critical end-points 
emerging from the zero-temperature plane (see FigHJ) . This bifurcation is accompanied by 
the appearance of a striking peak in resistivity |9| with curious, anisotropic dependence on 
the relative orientation of current, lattice and in-plane magnetic field When current 
flows in the crystallographic direction most parallel to the in-plane field, the resistivity peak 
rapidly decreases as the field is moved away from the c-axis. When it is nearly perpendicular 
to the in-plane field, the peak persists. Further indications of a "roof" delineating the region 
of anomalous phase behavior with field along the c-axis were provided by a kink in the 
longitudinal magnetization and a qualitative change in the temperature dependence of re- 
sistivity. Fig. CD uses new resistivity data to extend this roof in angle. Similar features occur 



elsewhere in the phase diagram 



111 ], with further bifurcations apparent upon approaching 



the ab-plane. These show a smaller resistance anomaly, but have the same characteristic 
anisotropy. 

Beginning with a heuristic discussion of the physics of the LOFF state and its magnetic 
analogue, in the following, we will describe how the the Wohlfarth-Rhodes jl^] band picture 
of metamagnetism is extended to allow the possibility of spatially modulated magnetic 
phases. In order to deduce the effects upon the broader phase diagram, we turn to a 
Ginzburg-Landau expansion of the microscopic Hamiltonian. The key physics is revealed in 
an expansion along the line of metamagnetic critical end-points through a vanishing stiffness 
to spatial modulation of the transverse magnetization. This leads to a reconstruction of the 
phase diagram. Finally, we describe how our picture explains the behaviour of SrsRu207 — 
capturing both the experimental phase diagram and the properties of the anomalous phase. 

A BCS superconductor is formed by binding electrons at the Fermi surface with opposite 
spin and momentum (k, f and — k, j) to form Cooper pairs. A magnetic field imposes 
a Zeeman energy cost on the superconductor which is balanced against the condensation 
energy. When Zeeman energy dominates, the superconducting state is destroyed; Cooper 
pairs are broken allowing a spin polarization to develop. The transition from a superfluid 
to a normal phase can occur via an intermediate inhomogeneous condensate, the LOFF 
phase 

HQ. 

By pairing electrons into a state with non-zero total momentum (k + q/2, f 
and — k + q/2, j), the reduction in condensation energy due to modulation is offset by a gain 
in Zeeman energy. The precise texture of the superconducting order depends sensitively 
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upon microscopic details j^]. 

A similar mechanism can apply to itinerant magnets: A spatially-modulated magnetic 
phase may intervene between the high- and low- magnetization states of a metamagnet. To 
form a ferromagnet, there must be an energetic gain in transferring an electron from a spin- 
down to a spin-up state of the same momentum. In a Stoner model, this is due to Coulomb 
exchange energy acquired at the expense of kinetic energy. Extending the Stoner model 
to include a band dispersion with peaks in the electronic density of states (DoS) leads to 
metamagnetism [ij], As the Fermi surface of, say, majority carriers approaches its van 
Hove filling, the single-particle energy cost in changing its filling is reduced. This can lead 
to a step change in the magnetization at certain values of the external field. 

Inhomogeneous mag netic states can be stabilized by peaks in the DoS in a similar way to 
spin density waves [lT], Q]- The simplest inhomogeneous phase formed from a ferromagnet 



is a spin spiral 19]. A spiral of the right wavevector distorts the Fermi surface so that some 
regions are brought closer to their van Hove filling (see Fig. [2]). The reduction in single 
particle energy cost due to occupying states near to the peak in the DoS can outweigh the 
single-particle energy costs from elsewhere. This leads to peaks in the transverse magnetic 



susceptibility 17| and ultimately provides a mechanism by which a metamagnetic transi- 
tion can split the transition between low and high magnetization occurring via a phase of 
inhomogeneous transverse magnetization. 

Such behavior can be shown explicitly for a Stoner model with band dispersion: TC = 
Y^k a =i i e k^k,o- — Un^hi —fiBH(n^ —ni), where is the number operator for electrons with 
momentum k and spin a and is the total number operator for spin-up and spin-down 
electrons, respectively, is the electronic dispersion — we focus, without loss of generality, 
upon a two-dimensional tight-binding dispersion with next-nearest neighbor hopping. As 
noted above, this model displays metamagnetism [ni [li]]. 

Inhomogeneous phase formation leads to a reconstruction of the metamagnetic phase 
diagram that is best revealed through a Ginzburg-Landau expansion. The thermodynamic 



properties of a metamagnetic system can be developed as a Landau expansion in magneti- 



zation density, M, as 
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f3F L = rM 2 + uM 4 + vM 6 - h ■ M , (1) 
where h = heu is the external magnetic field. The coefficients of this expansion for the 
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FIG. 2: a) Energy contours for a next-nearest neighbor tight-binding model with Fermi surfaces 
for minority and majority electrons shown in blue and red, respectively (left). The minority and 
majority bands are distorted by a spiral modulation with non-zero transverse magnetization and 
wavevector (right), b) DoS of minority and majority states with a uniform magnetization (left) 
and with a spiral distortion (right). In the former, the Fermi surface lies just below a peak in the 
DoS and in the latter it lies between two split peaks. 



Stoner model may be calculated in a standard manner following a perturbative expansion 
in interaction 2l|], and are specific functions of the external parameters; magnetic field, 
temperature and distance from van Hove filling. The point r = u = h = denotes the 
position of the parent tricritical point where the line of continuous Stoner transitions at zero 
field (r = h = 0, u > 0) bifurcates symmetrically into two lines of metamagnetic critical 
end-points parameterized by the conditions, 3mJ~l = 9 2 m Tl = d\i^L = 0. 

As we are interested in a reconstruction of the metamagnetic transition, it is convenient to 
shift our expansion from zero magnetization to an expansion about the magnetization along 
the line of metamagnetic critical end points. Setting M/M = (1 + <t>(r)) e\\ + </>j_(r), where 
M denotes the mean- field magnetization along the metamagnetic line 22j, and substituting 
into (JTl) gives 



hM 



(2) 



H and R parameterize deviations from the metamagnetic critical end-point perpendicular 
and parallel to the first order line. The dependence of the higher order coefficients on H 
and R can be neglected. Although we are interested in reconstructions of the metamagnetic 



transition that may lie outside the radius of convergence of the parent Landau theory (pQ), 
it turns out that an explicit expansion for the Stoner model gives the same coefficients as 
(J2J) when constrained to lie along the parent line of critical end points 23]. 

To allow for inhomogeneous phase formation, we consider a minimal gradient expansion 
of the free energy: 

(3F GL = (3F L + (K ± + K x 4> + K 2 <f) 2 + K 3 <f> 2 ± ) (V0 ± ) 2 

+L ± (V 2 </> ± ) 2 , (3) 

where the parameters K±, K 2 , K3 and L± are functions of the external parameters fixed by 
the microscopic theory. We have neglected gradient terms associated with cf). While such 
terms can lead to a spatial modulation, they do not lead to the phase reconstruction that 
we find. Gradient terms of fourth order and higher ought strictly to respect the lattice 
anisotropy 24j]. We consider the isotropic case for simplicity. 

The key ingredient introduced by explicit evaluation of the coefficients of the gradient 
expansion for the Stoner theory - that cannot be anticipated on purely symmetry grounds 
- is that K_l changes sign along the line of metamagnetic critical end points 25j. This 
indicates an instability to the formation of a spiral transverse magnetization. As thi s sp iral 
order is established, the effective 4 term changes sign leading to a tricritical point [10(] . 

Inhomogeneous phase 

Line of 
critical 
endpoints 




Metamagnetic sheet 

FIG. 3: Phase diagram for the Ginzburg-Landau theory with possible spin texture. Green sheets 
represent first-order transitions in eft. Blue sheets represent continuous transitions into the inho- 
mogeneous phase. The possible spin texture is constructed from four spin helices arranged in a 
square. The longitudinal magnetization has been supressed in this picture for emphasis. 



The resulting phase diagram is shown in Fig. [HJ The metamagnetic 
dislocated (symmetry broken) tricritical point as shown in green lOj 
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FIG. 4: The phase diagram rotated into the experimental orientation. The dashed line shows a 
trajectory through the inhomogeneous region. The inset shows the variation of longitudinal and 
transverse magnetization through this trajectory. 



embrace a region of inhomogeneous transverse magnetization in accord with our heuristic de- 
scription. This region is further enclosed by a surface of continuous phase transitions, shown 
in blue, at which the transverse magnetization falls to zero. The longitudinal magnetization 
shows a kink on this surface — a ghost of the transition in the transverse magnetization. 

The inhomogeneous magnetic structure may consist of a superposition of several wavevec- 
tors. The sum of these wavevectors must be zero to avoid a spontaneous spin current. A 
4-fold lattice symmetry (as in S^R^Oy) suggests four preferred wavevectors. There are 
two ways to superpose these: in pairs of ±q leading to a spin density wave in one of two 
directions that breaks the 4- fold rotational symmetry to 2; a superposition of all four sym- 
metry related wavevectors leading to a spin crystal which preserves the lattice symmetry. 
An example of the latter case is shown in the inset to Fig. [31 

Comparison with the experimental phase diagram, Fig. [TJ is obtained by expressing R, 
H and K± as functions of the experimental parameters T, 9 and h. These functions are 
expected to be analytic (as is confirmed by their detailed microscopic calculation) and in 
the usual spirit of the Ginzburg-Landau expansion their leading dependence near to the 
critical point is linear [261 ] . Here we choose to expand about the point along the line of 
metamagnetic critical end-points where K± changes sign. Fig. H] shows the result of such a 
correspondence. The natural parameters of our microscopic theory are field, temperature 
and band filling. An additional mechanism is required to map from filling to angle. One 
candidate is spin-orbit coupling [271] (which leads to an angle dependent Zeeman coupling) 
together with orbital effects of an in-plane field in a bilayer system. As the anomalous 
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behaviour of Sr 3 Ru 2 7 only appears in the cleanest samples, its origin must be sensitive to 
disorder. Our mechanism shows this sensitivity, since disorder smooths out features in the 
DoS. 

Spatially inhomogeneous magnetic structures lead inevitably to enhanced scattering in 
certain directions. In order to fully explain the anisotropy, there must be a mechanism for 
an in-plane magnetic field to align the magnetic inhomogeneity. Our simple model does 
not contain such a mechanism. We suggest that its origin lies in a modification to the 
dispersion due to in-plane magnetic field, which breaks the symmetry between different 
orientations of the underlying helices. In the anomalous phase this magnetic inhomogeneity 
leads to enhanced resistivity. With a magnetic field in the c-direction, the inhomogeneity 
does not break the crystal symmetry (at least macroscopically) and resistivity is isotropic. 
As the field is rotated into the plane, the magnetic inhomogeneity no longer preserves the 
lattice symmetry — either through the formation of an anisotropic spin crystal or by a 
preponderance of domains of spin density waves of one orientation. This anisotropy is 
reflected in resistivity. 

Spatial modulation of magnetization should show up as Bragg peaks in elastic neutron 
scattering in the anomalous region. Unfortunately, no such data exist. There are, however, 
pseudo-elastic data outside of the anomalous region consistent with fluctuations that would 



28]. 



freeze into the type of spin-crystals that we predict 

The mechanism of inhomogeneous magnetic phase formation presented here contrasts 
with two other proposals: i. Spin orbit interactions in systems without a centre of inversion 
symmetry lead to a Dzyalosinskii-Moriya interaction [29( that favors the formation of mag- 
netic spirals j^] and possibly magnetic crystals 31 1. We restrict attention to systems, such 
as Sr 3 Ru207, that have a centre of inversion symmetry, ii. Analysis of quantum fluctuation 
corrections to the theory of itinerant magnets suggests that they can induce metamagnetism 
and magnetic inhomogeneity 32j. Whether such effects are important in SrsRu207 is un- 
clear. We expect that van Hove singularities are characterized by larger energy scales and 
provide the dominant mechanism. Others have speculated that the anomalous phase in 



Sr 3 Ru207 may be a nematic metal with a d-wave distortion of the Fermi surface [33]. The 
topology of the phase diagram resulting from this distortion should be similar to ours if 
extended in angle. The main distinction is in the spatial modulation that we predict, which 
could be probed directly by neutron scattering. 



S 



In conclusion, it has long been established that the Stoner model with a peak in the 
DoS can lead to metamagnetism. We have shown that a generic band dispersion leads to 
a bifurcation of this metamagnetism by the intervention of a phase of spatially modulated 
magnetism analogous to the superconducting LOFF state. This behaviour might have been 
seen already in Sr3Ru207. Indeed, our analysis is rather general and its results may have 
broader applicability, e.g. NbFe2 34J exhibits a peak in resistivity associated with the 
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bifurcation of a metamagnetic transition and finite wavevector magnetic order and ZrZn 2 
may show similar features. 
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Lonzarich and A. P. Mackenzie for insightful discussions. 
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